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Abstract 

We derive exact relations that connect the universal C/fc 4 -decay of the momentum dis- 
tribution at large k with both thermodynamic properties and correlation functions of 
two-component Fermi gases in one dimension with contact interactions. The relations 
are analogous to those obtained by Tan in the three-dimensional case and are derived 
from an operator product expansion of the one- and two-particle density matrix. They 
extend earlier results by Olshanii and Dunjko [l[ for the bosonic Lieb-Liniger gas. As an 
application, we calculate the pair distribution function at short distances and the dimcn- 
sionless contact in the limit of infinite repulsion. The ground state energy approaches a 
universal constant in this limit, a behavior that also holds in the three-dimensional case. 
In both one and three dimensions, a Stoner instability to a saturated ferromagnet for 
repulsive fermions with zero range interactions is ruled out at any finite coupling. 
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1. Introduction 

The study of a non-relativistic system of fermions with spin-independent two-body 
interactions appears as a generic many-body problem in different areas of physics. Except 
for the particular case of one dimension (ID) 0, 0, [|, there are, unfortunately, very 
few exact results on this problem beyond the perturbative regime, which is typically 
not the one that is realized in nature. It is therefore of considerable interest to derive 
relations for the many-body problem that hold independent of the interaction strength. 
A well known example is Landau's Fermi liquid theory, which provides exact results 
for the low-energy excitations and the resulting response functions in terms of a few 
phenomenological parameters. The validity of these relations, however, requires that the 
interacting ground state is adiabatically connected to that of free fermions, an assumption 
that can hardly ever be justified starting from the microscopic Hamiltonian. Remarkably, 
considerable progress and insight into the fermionic many-body problem at arbitrary 
strength of the interactions has been achieved recently by Sinha Tan [5j in the particular 
case of zero range interactions. In this special case, it turns out that the momentum 
distribution exhibits a universal C/k decay as k approaches infinity. The constant C is 
the same for both particle species [H| and is called the contact, because it is a measure 
of the probability that two fermions with opposite spin are close together [B| (for a 
comprehensive recent review of the Tan relations see [6j). The contact determines the 
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change of the energy with respect to the interaction strength by a Hellman-Feynman like 
relation, the Tan adiabatic theorem Q. It also allows to calculate the energy from the 
momentum distribution Q. A crucial feature of the Tan relations is the fact that they 
apply to any state of the system, e.g. both to a Fermi liquid or a superfluid state, at zero 
or at finite temperature and also in a few-body situation. The only change is the value of 
the contact C. The origin of this universality was elucidated by Braaten and Platter [sj 
who have shown that the Tan relations are a consequence of operator identities that 
follow from a Wilson operator product expansion of the one-particle density matrix. 

Our aim in the present work is to derive the analog of the Tan relations for two- 
component Fermi systems in ID. At first sight, this appears to be of little interest because 
ID fermions with zero range interactions can be solved exactly by the Bethe-Ansatz 
[flLUflL 11 1 . As will be shown below, however, the Tan relations in ID provide information 
on observables that are not easily calculable from the Bethe-Ansatz, like the momentum 
or the pair distribution function. Moreover, they also apply at finite temperatures or in 
the presence of an external confining potential, where the Bethe-Ansatz fails. This is of 
particular relevance for the case of ultracold fermionic atoms, for which a model with 
contact interactions in fact provides a realistic description of atoms confined into a ID 
quantum wire geometry 12]. In particular, the use of Feshbach resonances in this context 
allows to tune the interaction strength simply by changing a magnetic field. This gives 
access to the whole regime from weak interactions to the limit of infinite attraction or 
repulsion. In the 3D case, the latter limit is reached in the unitary Fermi gas, where the 
two-body scattering length a diverges. The attractive branch of the unitary Fermi gas 
has a superfluid ground state and has been studied quite extensively in the context of the 
BCS-BEC crossover problem [l3|, [12, [Hj] • It provides an example of a non-relativistic 
field theory that is both scale and conformally invariant at kpa = oo 15, 1(|. The nature 



of the ground state on the - metastable - repulsive branch, in turn, is still unknown beyond 
the perturbative limit which - by construction - is a Fermi liquid [17| . In particular, 
it is an open question whether the repulsive 3D gas exhibits a Stoner type instab ility 
to a ferromagnetic state, as indicated by a renormalized Hartree-Fock calculation [18j. 
Experimental support for the existence of such an instability has been provided by the 
observation of a sudden decrease in the three-body loss rate and a minimum in the kinetic 
energy of the gas at a critical value kpa — 1.9 ± 0.1 of the interaction parameter [lij ]. 
Recent variational [2(| and numerical calculations 21, 22 1 in fact find that the ground 
state energy of the unpolarized Fermi liquid state exceeds that of a saturated ferromagnet 
for sufficiently strong repulsion. The associated value of the critical coupling constant 
turns out to be much smaller than that observed experimentally, a discrepancy that might 
be explained by the fact that the repulsive branch is only metastable. By contrast, as 
will be shown below, even in an equilibrium situation where the lifetime of the repulsive 
branch is assumed to be infinite, a combination of the Tan relations and a variational 
argument rules out a Stoner instability to a saturated ferromagnet for fermions with zero 
range interactions in both one and in three dimensions. 

Our derivation of the Tan relations is based on the Operator Product Expansion 
(OPE), which was developed independently by Wilson [23j and Kadanoff 124]. In fact, 
this method was implicitly used in this context first by Olshanii and Dunjko [lj in a study 
of the short distance behavior of the one-particle density matrix of the Lieb-Liniger gas, 
a system of ID bosons with zero range interactions [251 ] . Specifically, they found that 
the momentum distribution of this bosonic system decays like Cs/k A for large k. The 
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associated contact coefficient C'b turned out to be proportional to the pair distribution 
function g^ (0) at vanishing distance between the particles. As will be shown in section 
4 below, this is effectively the ID analog of the Tan relation for the asymptotic decay of 
the momentum distribution for bosons. It carries over to the two-component Fermi gas 
with only slight modifications. 

The paper is organized as follows: in section 2 we introduce our model and discuss 
its range of validity. The Tan adiabatic theorem is derived in section 3 as a simple 
application of the Hcllman-Fcynman theorem. In section 4 the OPE of the one-particle 
density matrix is used to connect the asymptotics of the momentum distribution with the 
pair distribution function for opposite spins at vanishing separation, which is essentially 
the contact in one dimension. Thermodynamic relations that connect energy and pressure 
with the contact are derived in section 5. An OPE of the density correlation function 
in section 6 shows that the contact also appears as a non-analytic contribution to the 
pair distribution function at short distance, giving rise to sum rules and a power law 
asymptote in the static structure factor. Explicit results for the contact at arbitrary 
interaction strengths are given in section 7 for the particular case of the balanced Fermi 
gas at zero temperature on the basis of the known Bcthc-Ansatz solution for the ground 
state energy. Finally, in section 8, we show that irrespective of an explicit solution of the 
many-body problem, a combination of the Tan relations and a variational argument rules 
out saturated fcrromagnctism for zero range interactions at any finite coupling, both in 
one and also in three dimensions. Details of the diagrammatic expansion for the OPE 
are presented in an Appendix. 



2. The ID Fermi gas with contact interactions 

We consider a Gaudin-Yang model for a two-component system of fermions in ID 
which interact with a two-body <5-function potential. Due to the Pauli principle, this 
interaction only affects fermions with opposite 'spin' a — J \ 1 ^. In the case of ultracold 
atoms, the two spin states might be two different hyperfine states of 6 Li or 40 K or they 
might describe a fermionic mixture with different masses m a . In a second quantized 
form, the Hamiltonian is expressed in terms of field operators ip a (R) which obey the 
usual anti-commutation relations. Including an external one-body potential V(R) and 
using units in which h = 1, the Hamiltonian H = J dRH(R) can be written in terms of 
a Hamiltonian density 

n ( R ) = J2 ^d R ^d R MR) + 9 ^W-^-MR) +Y;V*miMR)- (i) 

The (5-function interaction is a proper potential in ID and thus needs no renormalization. 
Its coupling constant g\ is conveniently expressed in terms of a one-dimensional scattering 
length ai by g\ = —l/m r a\, where m r is the reduced mass of the two interacting particles. 
In order to discuss the physical meaning of this scattering length and the question under 
which conditions a realistic two-body interaction potential can be replaced by an effective 
contact interaction, we consider the low-energy limit 

r(fc) = i + icot(<s(k)) ~ i + ik ai + o ((ke ± ) 3 ) (2) 
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Figure 1: Integral equation for the scattering of a pair of fermions. The big gray blob represents the 
amplitude iA{E). 

of the reflection amplitude that describes two-particle scattering in one dimension quite 
generally in a situation, where the particles are confined to the ground state of their 
transverse motion with size £± 12|. For a gas with total density n, the momenta for 
two-body scattering are typically of the order of the Fermi momentum kp = nn/2. The 
reflection amplitude can thus be replaced by its low-energy limit —1/(1 + ikai) provided 
that kp <C a±/£j_. In a regime of low densities, therefore, any two-body interaction for 
which the low-energy expansion in equation ^) holds, can be replaced by an effective 
contact interaction giS(x) with coupling constant g\ = —\/m r ai, for which the reflection 
amplitude equals —1/(1 + ikai) at arbitrary values of k. Specifically, for ultracold atoms 
whose interaction is described b y a 3D pseudopotential with scattering length a, the 
resulting value of a\ is given by [261 ] 



ai(a) = — ± + A£ x , (3) 



where A = — £(l/2)/-\/2 ~ 1.0326 is a numerical constant. Note, that a negative value 
ai < of the ID scattering length corresponds to repulsive interactions and vice versa, a 
situation that is precisely opposite to the 3D case. For a± > 0, the scattering amplitude 
has a pole at k = ija\ with a positive imaginary part. This describes the two-body 
bound state of an attractive (5-function potential, with a\ the size of this bound state. 

For the calculation of the diagrams that appear in the Wilson-OPE for the one- and 
two-body density matrix in sections 4 and 6 below, it is convenient to use the amplitude 
A(E), which is related to the on-shell T-Matrix at the center of mass energy E = k 2 /2m r 
by A(E) = —T(E). It obeys the Lippmann-Schwinger equation shown diagrammatically 
in figure [TJ For the 5-function interaction, the Lippmann-Schwinger equation is solved 
analytically by 

AE) = ^-r(k) = ±- 1 i (4) 

m r TTl r a\ . 

which is just the associated refection amplitude up to a factor —ik/m r . 

3. The Tan adiabatic Theorem 

For ID systems with zero range interactions, the expectation value of the interaction 
part of the Hamiltonian 

(H')= 9l [ dR(^l^(R)) (5) 



is quite generally determined by the value of the pair distribution function for opposite 
spins 

5 g ) ( J R,i?') = (^(i?)^(i?')^(^)^(^))/K(^)«i(^)) (6) 
at vanishing distance R = R' between the particles, as noted by Lieb and Liniger in the 
Bose gas case (25j . The interaction energy in ([3]) is therefore simply g\ times an integral 
of the probability density for two fermions with opposite spin to be at the same point in 
space, a quantity that is proportional to the contact density discussed in section 4 below. 
More precisely, it is convenient to define the contact C as an extensive variable by 

C = 4 gi m 2 r (H') = Aglm 2 r J dR (^\^{R)). (7) 

This variable will play a central role in the following. It depends both on the microscopic 
parameters like the scattering length a\ and also on the specific state under considera- 
tion. While the contact is basically the expectation value of the interaction part of the 
Hamiltonian, it also determines the variation of the total energy E with the scattering 
length oi. Indeed, by the Hellman-Feynman theorem, we have 



dE _ /dH_\ 



1 1 



T~ = ( TT ) = -— / dR WWWMR)) ■ ( 8 ) 
da\ \ oai I a\ m r J + 

Using the definition of the contact, this immediately gives the ID analog of the Tan 
adiabatic theorem Q in the form 

±E = %2l. (9) 
dai 4m r 

A knowledge of the dependence of the contact C on the scattering length therefore 
determines the associated rate of change of the total energy. Note that the derivative in 
© has to be taken with all other variables constant, in particular the entropy if thermal 
states are considered. This will play a crucial role in the derivation of the pressure 
relation in section 5. A particularly simple form of the Tan adiabatic theorem is obtained 
by considering a translation invariant situation with vanishing external potential V = 0. 
Applying the Hellman-Feynman theorem to the Hamiltonian density in equation ([T]) 
then shows that the pair distribution function glJ(0) for opposite spins at zero distance 
is determined by the derivative of the energy per length with respect to the coupling 
constant g\ via 

^= W f>). (10) 

This relation will be used in section 7 below to determine the contact explicitly in the 
ground state of the balanced two-component Fermi gas from the known Bethe-Ansatz 
result for the ground state energy for arbitrary coupling strength. 



4. Asymptotics of the Momentum Distribution 

In the following, we want to show that the contact introduced in (J7J determines not 
only the interaction energy and the derivative of the total energy with respect to a\ but 
also appears in the asymptotic behavior of the momentum distribution 

n a (k) = (^(fc)Vv(fc)), (11) 
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where ip a (k) is the Fourier transform of the field ip a . The momentum distribution n a {k) 
is normalized such that the integral over all momenta gives the total number N a of 
particles of type a. This will be more convenient for the following arguments compared 
to the standard intensive normalization of n a (k) to the respective densities n a . In terms 
of the fields tp^(R) and ip cr (R), the momentum distribution can be expressed by 



Its behavior at large momenta k can therefore be obtained from an expansion of the 
one-particle density matrix (ip^(R)ip a (R + x)) for small separations x — > 0. To do this in 
practice, we use the Wilson-Kadanoff Operator Product Expansion (OPE) which states 
that the product of two local operators separated by a short distance, can be expanded 
as the sum of local operators with coefficients that depend only on the separation. For 
the specific case of the one-particle density matrix, this statement reads 



where the local operators O atn (R) are products of the quantum fields and their deriva- 
tives. The Wilson coefficients c atn (x) encode all the short distance behavior, in particular 
possible singularities. The expansion is to be understood as an asymptotic one, i.e. for 
any positive integer / only a finite number of terms on the right hand side vanish more 
slowly than \x\ l as x — > 0. Knowledge of the leading terms in the OPE determines the 
large k behavior of the momentum distribution via 



Since (|13[) is an operator equation, it holds for matrix elements between arbitrary 
states. This allows to determine the Wilson coefficients ^^(x) by a matching procedure. 
We compute the matrix elements of the left hand side and right hand side of equation 
dJ3]) between states (x| , |</>), expand the left hand side with respect to x, and demand 
both sides of the equation to be equal up to some order in x. The Wilson coefficient of 
an operator O a ^ n {R) may be determined from the matrix element between any states, 
for which the matrix element (\\ Ca,n{R) |</>) is non-zero. Except for this constraint, one 
may choose the simplest possible states one can think of. 

Our explicit calculations follow closely the method used by Braaten and Platter in the 
three-dimensional case 8,27]. For the coefficients of one-particle operators such as iplipa , 
we will choose one-particle plane wave states (p' | , \p) containing one particle of the species 
a with momentum p' or p respectively. For two-particle operators, such as ^jV'JV^V'T' we 
choose two-particle scattering states (±p'|, |±p) in the center of mass frame, containing 
one particle of each species with momenta ±p and ±p' respectively. The computation of 
these matrix elements will be simplified by drawing the corresponding diagrams and using 
the Feynman rules for the theory and the operator vertices given in |Appcndix A| We 
will consider the OPE of i/j^.(R)ijj^(R + x ) to be able to draw the diagrams in a definite 
way - the matching for the one-particle density matrix of species a =\. is completely 
analogous. The matrix elements for the matching of the one-particle operators can be 
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Figure 2: Diagrams for the matrix elements of one-particle operators between one-particle states. 
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Figure 3: Diagrammatic contributions to the left hand side of the OPE of the one-particle density matrix, 
using two-particle scattering states. 



represented by the diagrams shown in figure [3J In particular, the matrix element of the 
one-body density matrix in single-particle plane wave states is equal to 



(p'\4(R)iJ t (R + x) \p) 



-p')R e ipx 



(15) 



The right hand side consists of operators of the type ^(du) m if!^, whose matrix elements 



are given by 



(p'\4(d^)(R) \p) = (i^e'f"'' 8 



(16) 



Comparing both sides of the OPE, the Wilson coefficients of the local one-particle oper- 
ators 7/>j(<9R) m -0t have to match the exponential term e lpx and are thus given by x m /m\. 
The contribution coming from these operators is analytic in x and is just equivalent to a 
Taylor expansion of the one-particle density matrix around x = 0. Note the slight differ- 
ence compared to the matching performed in [§] , because we have put the x-dependence 
only in the argument of the i/v-field. Thus, there are no analytic contributions involving 
derivatives of in our case. In analogy to the three-dimensional case there are, however, 
also non-analytic contributions in the short distance expansion of the one-particle density 
matrix. The leadingterm is, in fact, again connected with the operator i/jtip]ip^ip-f that 
also appears in 3D [8|. To show this, we calculate the matrix element of the one-particle 
density matrix between the two-particle scattering states. Diagrammatically, this can be 
represented by the sum of the four diagrams shown in figure [3) The first three diagrams 
with one or no scatterings contribute purely analytic terms. Together with the analytic 
part of diagram |3(d) | they are matched by the one-particle operators whose Wilson co- 
efficients are already known. The high momentum region in the loop integral gives rise 
to a non-analytic contribution to the diagram [3(d)] of the form 



im 2 r A(E)A(E')- 



1 



\P'\ 2 
7 



i\p\\x\ 



(17) 
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Figure 4: Diagrams for the matrix element of the operator i/>|i/>|i/i^i/>^ 

where A(E) is the scattering amplitude given in equation (Q} while E = p 2 /2m r and 
E' = p 2 /2m r are the associated energies. The expansion of the bracket in equation 
(|17[) in powers of \x\ has a leading order non-analytic term proportional to \x\ 3 , which, 
after Fourier transform, will give rise to a l/k 4 tail in the momentum distribution. 
To determine the Wilson coefficient of the operator ipj^ipjjip^ipf, we compute its matrix 
element between the two-particle scattering states, which is given by the sum of the four 
diagrams shown in figure 3] Using the Lippmann-Schwinger equation in the form 

-m = 1 + iA(E) ^ ( 18 ) 
5i \P\ 

where the factor m r /\p\ is the value of the loop integral, the sum of the four diagrams 
can be reexpressed and the matrix element is 

(±p'\ i/tfyfyMR) \± P ) = A ^f E '\ (19) 

9i 

This matches the |a;| 3 -term in the expansion of (fTF)) with Wilson coefficient g 2 m 2 .\x\ 3 /3. 
In summary, the short distance OPE (|13[) of the one-particle density matrix, up to 
order x , is the sum of three analytic contributions and a leading order non-analytic 
contribution proportional to \x\ 3 . The operator responsible for this non-analyticity is 
ij^ip^ip±ipf(R) , in complete analogy to the three-dimensional case Q. Explicitly, the 
Wilson-OPE for the one-particle density matrix up to order x 3 reads 

V£(fl)^(i2 + z) = X;^-^(W (20) 

m— 

Clearly, the zeroth order term at x — is just the density of species a at location 
R, as expected. The second order Wilson coefficient is associated with the operator 
0<t,2(-R) = V4- \dji4>o-) (R) which is the kinetic energy density. Expanding the one-particle 
density matrix in the ground state of a uniform, balanced Fermi gas as (ipl-(x)ip a (0)) = 
nil + C2n 2 x 2 + ...)/2, this contribution can be expressed in terms of the dimensionless 
ground state energy e("j) introduced in (|52p below in the form 



C2 = -[e(7)-7e'(7)]/2, 
8 



(21) 



in close analogy to a result derived by Olshanii and Dunjko [l[ for the Lieb-Liniger gas. 
Since this is an analytic contribution in x, it does not give rise to a power law in the 
associated momentum distribution. 

The leading non-analytic contribution appears only at third order in x unlike in three 
dimensions, where it is of order r = |x|. The absence in ID of a contribution proportional 
to |x| may be traced back to the fact that both the kinetic and the interaction energy 
are finite in ID. A term ~ \x\ in the one-particle density matrix would lead to a l/fc 2 -tail 
in the momentum distribution, thus making the kinetic energy divergent. The leading 
term in the momentum distribution for high k can be obtained by inserting the OPE 
([2"U|) up to order |x| 3 into equation (TH)) . The resulting power law tail is obtained from 
the formal Fourier transform 



dx\x\ a e~ ikx = 2T{ a +l)co S -(a + 1) (22) 



7T 



1 



\k\ a + i: 



of \x\ a , with a a non-even real number. Specifically, for a = 3 one obtains 12 /k 4 . With 
the prefactors given by (|20[) , we have thus shown that the momentum distribution 

M*) ^ § (23) 

of a two-component Fermi gas in one dimension with contact interactions, exhibits a 
power law decay C/k A at large momenta k. (In practice, for a uniform gas with ID 
density n, the asymptotic behavior shows up at k 3> n). The constant prefactor C, called 
the contact, is the same for both species and is positive by definition. Microscopically it 
is given by the expression Q which is proportional to the expectation value (H') of the 
interaction energy. In the translational invariant case, (H') ~ L is extensive since the 
pair distribution function ((5]) only depends on the coordinate difference. The resulting 
intensive contact density C — C/L is then given by 

C = A ntn ^(2)( ). (24) 



This result is completely analogous to the expression 



C e = 4«V 2) (0) (25) 



for the contact density of the ground state of the uniform Lieb-Liniger gas derived by 
Olshanii and Dunjko [lj. Note that in the non- interacting case gi = 0, the non- analytic 
terms disappear and there is no tail, since the momentum distribution is just a Fermi- 
Dirac or Bose-Einstein distribution for free particles which exhibits no 1/fc 4 tail at any 
temperature. 



5. Energy and Pressure Relation 

In the following we will show that the connection between the asymptotic behavior 
of the momentum distribution and the interaction part of the Hamiltonian also allows 
to derive three other universal relations for the two-component Fermi gas with contact 
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interactions which are again analogous to those in the 3D case. We start with the energy 
relation: 

The total energy E = (H) = (T) + (I) + (V) of the Fermi gas in the presence of an 
additional one-body confining potential V(R) can be expressed in terms of the contact 
C and the momentum distribution n CT (fc), as 

E = S/|^M*)-» I ji ; + (V). (26) 

(T 

This relation follows quite simply by inserting the definition ([7]) of the contact into the 
Hamiltonian (JXJ) , and rewriting the kinetic energy in momentum space 

dR,{d R iPl(R)d R MR)) = I Y e ! ($(*Mt(*0) • (27) 

^ V ' 

n„(fc) 

The energy relation expresses the energy as a functional of the momentum and - in the 
inhomogeneous case - also the density distribution. It serves as a basis from which, 
together with the adiabatic relation, other Tan relations will follow. A special case of 
the energy relation is the virial theorem: 

For an harmonic trapping potential V(R) — mui 2 R 2 /2, the energy can be written as 
the sum of the trapping energy and the contact in the following way: 

E = 2 { V)- ai ^-. (28) 

This relation can be derived by a simple scaling argument. The only three energy scales 
are lu, l/m r a\ and ksT. The free energy 

F(T, N t , Ni,w, at) = uf (—, zr^s,^ N x ) (29) 

can therefore be expressed in terms of a dimensionless function /, which depends on the 
ratios fcsT/w, w/(l/ma 2 ) and the two particle numbers N^,N±. From (|29p. one can 
deduce the simple scaling law 

F{Xr,N^N i ,>^,\-^a l ) = \F{T,N t ,N i ,u,ax), (30) 

The derivative of equation (|30D with respect to A at A = 1 yields 

T m +u £;-\ a ^) F = F > (31) 

where all the partial derivatives arc to be understood as leaving all other system variables 
constant. Since the free energy is just the Legendre transform of the energy, its partial 
derivatives at constant temperature T with respect to uj and a\ are equal to those of the 
energy at the associated value of the entropy. Therefore, using OF/dT = —S, the energy 
turns out to obey the differential equation 

10 



This leads immediately to the relation (|28| by using the Hellmann-Feynman theorem 
with the trapping frequency w as a parameter of the Hamiltonian and the adiabatic 
relation ©. Again, the virial theorem is analogous to a corresponding relation in 3D that 
was derived first by Thomas et al. 28] for the unitary gas at a = oo and was extended 
to finite scattering lengths by Tan [29]. Note that at unitarity or - correspondingly - 
at a\ = in ID, the virial theorem allows to extract the total energy from the density 
profile n(R) in the harmonic trap. 

Finally, for the homogeneous (i.e. V a (R) = 0) two-component Fermi gas with contact 
interactions, there is a pressure relation which connects pressure and energy density £ 
(i.e. energy per length) by 

+ (33) 

where C = C/L is the contact density. For its derivation, we consider the free energy 
density T of the homogenous system. Since ep-f ~ n| and n^ep^ are characteristic scales 
for energy and energy density, dimensional analysis requires that T can be expressed in 
terms of a dimensionless function 

7"(T, ai, rif, nj,) = an^f (a T T/n^, n T ai, n ; /n T ) , (34) 

where n a is the density of species and > without restriction of generality. The 
function / depends on three dimensionless arguments which set the scale for temperature, 
interaction strength and a possible imbalance of the densities. The constants a and ap 
carry factors that are needed to make the arguments dimensionless and are irrelevant 
for the following. Apparently, equation (f34|) implies the following scaling-behavior of the 
free energy density 

T(X 2 T, A _1 ai, Xn-f, An;) = A 3 J"(T, at, rif, ni), (35) 

for an arbitrary, dimensionless parameter A. Taking the derivative of equation (|35p with 
respect to the parameter A at A = 1 yields 



d d d d 

oT on^ an^ oa\ 



T = 3J". (36) 



Now, quite generally, the grand canonical potential per length for a two-component 
system is given by 

J = -P = F - rif/Mf - . (37) 

This allows to replace the free energy density J- in equation (|36p by n^/i^+n^/i;— p. Using 
the thermodynamic relations dp J- = — s, where s is the entropy density, and d na J- = \i a 
results, after rearranging the terms, in 

d 

P = 2 (nfA*t + HHi -p + Ts) +oi ^— T. (38) 

v „ ' 0(1 l 

=T+Ps=£ 



Since the partial derivatives in equations (f36|) and (f38f are to be understood as leaving 
the other system variables constant, we can replace the d ai J- term by the adiabatic 
relation for the energy density. The resulting equation is the pressure relation ([3"3"]l . 
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Note that the sign of the correction term is opposite to that in the three-dimensional 
pressure relation (29[. At first sight, it therefore seems that a repulsive interaction a\ < 
lowers the pressure and vice versa. However, upon expressing 2£ in terms of a kinetic 
and an interaction energy part using the energy relation (|26p , we see that the interaction 
energy per length is equal to — aiC/4m r , which is positive for repulsive interactions, and 
negative for attractive ones, as expected. 

A quite nontrivial consequence of the pressure relation arises by considering the limit 
Oi = i.e. infinite interaction strength. In this limit equation Q33p predicts that 
pressure and energy density are related in precisely the same manner as for an ideal 
gas of non-relativistic particles in one dimension (note that this requires the product 
aiC(ai) ~ gl, (0)/ai to vanish in the limit a\ — > 0~, which is indeed the case as will 
be discussed in section 7 below). This surprising conclusion is a consequence of the 
fact that a ID Fermi gas with contact interactions is scale invariant at infinite interac- 
tion strength ai = 0. The underlying reason why an ideal gas type equation of state 
appears in this limit can be understood in simple terms by noting that the balanced 
two-component Fermi gas at infinite repulsion behaves like a free Fermi gas with a dou- 
bled value of the Fermi wave vector. Effectively, an infinite zero range repulsion in the 
two-component system is equivalent to a Pauli exclusion principle also between fermions 
of opposite spin. The resulting equation of state is thus p — 2£ at arbitrary temperatures 
despite the fact that it is a strongly interacting system. The same holds true in the case 
of infinite attraction ai — > + because the two-component Fermi gas in this limit is a 
Tonks-Girardeau gas of bosonic dimers, i.e. effectively an ideal single-component Fermi 
gas with kp — im/2 [3(J [Hj]. Note that in this case, both the energy density £ and the 
contact density C contain an infinite constant arising from the two-body binding energy, 
which diverges as g\ — > — oo. This constant, however, cancels in the equation of state, 
which at infinite attraction - is again the same as if there were no interactions at all. 



6. Short distance expansion of the pair distribution function 

In our discussion so far, we have only used the behavior of the one-particle density 
matrix at short distances. Its leading non-analytic contribution defines the contact and 
gives rise to Tan relations that are closely related to those in the 3D case [B|, 0, Q 
and also to those for ID Bose gases [![. In the following we will show that the contact 
also arises as a non-analytic contribution in the OPE of the two-particle density matrix 
which gives additional relations for the pair distribution function and the related static 
structure factor. Specifically we will derive the following short distance expansion 

Ar(*-fH(*+f)=(i- 

x 
+ 2 

of the f — i density-density correlator up to linear order in the separation x. The fact 
that the x — > limit of the pair correlator is simply the contact operator C(R) is not 
surprising. Note, however, that in three dimensions the relation is more complicated and 



2\x 



(39) 

n^d R fii - (dR-h^n^} (R) + 0(x 2 ) 
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Figure 5: Diagrams for the matching of the density— density correlator on the left hand side of its OPE. 
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C(R) 



ift^w (40) 

due to anomalous scaling of the contact operator. Equation (|39p shows that in ID, the 
contact is also responsible for the leading order non-analytic contribution to the pair 
distribution function, which is proportional to \x\. 

To prove the OPE given in equation (|3"5)l . first note that one cannot simplify the 
matching of some operators in one-particle states, as was the case for the OPE of the 
one-particle density matrix. Thus, we use scattering states (±p'| and |±p) to compute 
the matrix elements of the operators on both sides. The left hand side is the sum of the 
four diagrams shown in figure [5l Evaluating these yields 



i|p'IM 



-*(p~p')x +mrl ME')e- ipx ——+m r iA(E)e ip ' 

If] 



»»IpIM 



\P\ 



-+m 2 r iA(E)iA(E') 



e *|p|M g*|p'IM 



(41) 

Expanding this expression to linear order x, there are analytic contributions proportional 
to x° and x 1 plus a non-analytic one of order \x\. The cc°-term can be rewritten, using 
the Lippmann-Schwinger equation (fl8|) as A{E)A{E') / g\. The non-analytic |x|-term is 
given by 

2m r \x\ 



iA{E) [l + iA(E')j^- 



iA(E') l + iA(E)-^- 
\P\ 



where we have again used the Lippmann-Schwinger equation ([T 
and the |x| parts now reads 



(1 + 2m r g 1 \x\) 



A(E)A{E') 



A(E)A(E'), 
9i 

(42) 

The sum of the x° 



(43) 



Recall, that A(E)A(E')/gf is just the matrix element of the contact operator ^■0|'0^'0t C^) 
in these states, which is given in equation The Wilson coefficient of the contact op- 
erator, up to order \x\, is therefore given by (1 + 2m r gi\x\). We abbreviate the matching 
of the analytic part by noting that the matrix element of the operator h^dnh^ — {dfin^)n^ 
between the scattering states, which can again be represented by the diagrams in figure 
[U is given by 

2 -i(p-p') + (- ip )iA(E')j£- + ( tp ') t A(E)^ 
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(44) 



This matches exactly the linear term in the expansion of equation (l41~|) with Wilson 
coefficient x/2, thus completing the derivation of the OPE ((39|) of the density-density 
correlator. 

In order to discuss the consequences of the short distance singularity in equation (|39[) 
for physical observables, we consider the full pair distribution function g^'{R, R'). Quite 
generally, it is defined by the density-density correlation via 

(n(R)h(R')) = n(R)n{R')g {2) (R, R') + 6(R - R')n(R) . (45) 

Considering, for simplicity, the translation invariant case where the pair distribution 
function only depends on the separation x = R' — R of the two particles, the total pair 
distribution function 

■9 (2) (*) = ^ («t 9n (*) + n l 4? ( x ) + 2n t% 9n (*) ) (46) 

(21 

of a two-component Fermi gas splits into separate spin contributions g a a , (x) that all 
approach unity as x — > oo. For small separations x — > 0, the pair distribution function 
for equal spins will vanish like x 2 because the antisymmetry due to the Pauli principle 
forces the many-body wave- function to vanish linearly ~ (x, — Xj) as two coordinates 
Xi and x.j of fermions with equal spin approach each other. As a result, equation (|39[) 
implies that the total pair distribution function at short distances is given by 

g^x) = 2 ^g^0)(l-M + O { x^ (47) 

up to linear order in x. For repulsive interactions a\ < 0, therefore, the pair distribution 
function rises linearly from its value at zero separation. The associated short distance 
singularity ~ |x| gives rise to a 1/q 2 power law tail in the associated static structure 
factor 

S(q) = 1 + n [ dxe- iqx \g {2) (x) - ll . (48) 



Introducing a dimensionless coupling strength 7 = — 2/(na\), the Fourier transform of 
the non-analytic contribution in (|47[) leads to the asymptotic behavior in the form 

S(q^ 00) = l-7<4?(°) + (49) 

The tail of the static structure factor is therefore proportional to the contact density in 
equation (f24|) . A similar connection holds in 3D as pointed out by Hu et. al. [32]. In the 
3D case, however, the tail is proportional to C/q, which is a result of the l/r 2 -singularity 
at short distances described by equation (|4"0)) . Note that the tail in vanishes in the 
infinite repulsion limit 7 —J- 00 because g^, (0) ~ I/7 2 , as will be discussed in section 
7 below. At 7 = 00, the pair distribution function will in fact be equal to that of a 
non- interacting single-component Fermi gas with Fermi wave vector kp = 2kp = 7m, 
while the momentum distribution n a (k) will still have the nontrivial form of a Luttinger 
liquid (see section 7 below). 
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Finally, we note that in ID the contact density also arises in a sum rule for the static 
structure factor. Indeed, the Fourier inversion 

9 {2 \x) = l + - f p-e^[S(q)-l] (50) 
n J 2n 

of equation (|48p at x — immediately implies a sum rule 

*3%» = i + iff m-n <»> 

which connects the integrated static structure factor with the contact density. Here, we 
have used the fact that the contributions to the pair distribution function for equal spins 
vanish at x — and that S(q) is an even function. In contrast to 3D, where <7^(0) is 
not defined because even the zeroth moment of S(q) — 1 diverges due to the 1/g-tail, the 
finiteness of the pair distribution function at vanishing distance in ID guarantees that 
the deviation of the static structure factor from its trivial limit (q) = 1 that appears 
for both an ideal classical gas or an ideal Bose-Einstein condensate, is integrable. 



7. Contact density for uniform gases 

As mentioned in the introduction, the Tan relations that have been derived above are 
valid for an arbitrary state, both at zero and at finite temperature and in the presence 
of a non- vanishing one-body potential V(R). The simplest case of interest, for which 
explicit results for the contact density at arbitrary coupling strengths can be given on 
the basis of the Bethe Ansatz solution, is the ground state of a uniform gas with total 
density n, in which the two spin components are equally populated. In the repulsive case, 
this gas is known to be a Luttinger liquid [33l |J] . Quite generally, the one-particle density 
matrix of a Luttinger liquid of fermions decays with a power law ~ x~( a+1 ^ expifc^x at 
zero temperature. The resulting momentum distribution therefore exhibits a singularity 
~ \k — kp\ a at the bare Fermi wave vector kp, replacing the simple jump from one to 
zero of an ideal Fermi gas. The exponent a = (K — 1) /2K is related the Luttinger 
parameter K which - for the case of contact interactions - continuously varies between 
K = 1 for the non-interacting gas and K = 1/2 for infinite repulsion |4|. For attractive 
interactions, the balanced two-component Fermi gas is a Luther- Emery liquid 34, 30j . It 
has a finite gap for spin excitations due to the appearance of bound pairs with opposite 
spin, thus eliminating the singularity of the momentum distribution at hp. In both cases, 
the strength of the correlations in the ground state is characterized by a dimensionless 
coupling constant 7 = — 2/noi, which is inversely proportional to the density. In ID, the 
strong coupling limit \y\ 3> 1 is therefore reached at low densities. This somewhat coun- 
terintuitive fact can be understood by noting that low densities imply small momenta 
and, moreover, the ID scattering amplitude ([2]) has its maximum phase shift 5(0) = it/2 
as k — >• because ID potentials become impenetrable at zero energy, i.e. r(k —> 0) = — 1. 

The ground state energy of the uniform gas can be determined from the solution of 
the Gaudin-Yang integral equations, which can in fact be solved for arbitrary polarization 
(risf — ni) 7^ 35]. In the case of a balanced gas with equal masses, the ground state 
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Figure 6: The pair distribution function (0) at zero distance (black curve) and the associated dimen- 
sionless contact density s (orange curve) in the ground state of a balanced homogeneous Fermi gas as 
a function of the dimensionless coupling 7 = —2/na\. The black dashed line —7 on the left originates 
from the binding energy of two atoms, the right dashed curve follows from the asymptotic behavior of 
the energy on the repulsive side for large 7, given in 35]. The gray dashed line on the right is the 
asymptotic value s(oo) = 0.749. . . of the dimensionless contact density at infinite repulsion. 



energy per particle 

E h 2 n 2 

~N = ~~2m~ ' 6(7) (52) 

is conveniently expressed in terms of a dimensionless function e(7). For a repulsive gas 
with 7 > 0, the function e(j) increases monotonically from the value e(Q) = 7r 2 /12 that 
corresponds to the non-interacting two-component Fermi gas in ID to e(oo) = 7r 2 /3. In 
the attractive case 7 < 0, it is convenient to subtract a contribution — 7 2 /4 from the two- 
body binding energy. The relevant many-body energy §(7 < 0) = e(7 <_0) + 7 2 /4 then 
decreases monotonically from e(0) = 7r 2 /12 to e(— 00) = 7r 2 /48 [36|, [3(J |31| . In the case 
of a balanced gas, the connection (flO]) between the derivative of the energy per length 
with respect to the coupling constant g\ and the pair distribution function at vanishing 
separation reduces to the simple relation 

4f(0)=2e'( 7 ), (53) 

where e'{^) denotes the derivative of the dimensionless ground state energy defined in 
equation ((52]) . 

In figure |51 we show the pair distribution function as a function of the dimensionless 
coupling 7 for both the attractive and the repulsive gas as determined from a numerical 
solution of the Gaudin-Yang integral equations. Despite superficial appearance, the 
function is non-analytic at 7 = because the presence of pairing on the attractive side 
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7 < leads to an expansion of the form (0(—j) is the usual theta-function) 



4f(O,7<O) = l + | 7 |0(-7)-^ln 2 (| 7 |) + ... ■ (54) 

The contribution | 7 |0(— 7) arises from the two-body binding energy and dominates the 
behavior in the limit | 7 | 3> 1. The logarithmic term, in turn, is a non-analytic cor- 
rection to the mean field energy in the presence of pairing, whose many-body spin-gap 



A ~ exp — 7r 2 /2|7| is exponentially small as 7 — > 0~ [36|, |30(. On the repulsive side, the 
probability for two fermions with opposite spin to be at the same point in space continu- 
ously decreases to zero with increasing strength of the repulsion. In the limit 7 ^> 1 the 
decay is ~ 1/7 2 - As we will see, this reflects the fact that the contact density approaches 
a finite, universal value for a Fermi gas with infinite repulsion. To discuss the tail in the 
momentum distribution, it is convenient to define a dimensionless measure s of the value 
of the contact density C = s ■ kp by factoring out the Fermi wave vector kp = wn/2 as the 
characteristic scale on which the (intensive) momentum distribution n a {k) varies. For 
the ground state of the balanced gas, equation (|24l) then implies the simple connection 

s =±72. fl W( ) = A 7 V( 7 ), (55) 

between the dimensionless contact density s and the derivative of the ground state energy 
e(7). It is shown quantitatively in figure [6j again for both positive and negative values 
of 7. For a weakly repulsive (or also attractive) gas, s( 7 — !• 0) = A^ 2 /ir A vanishes as it 
should, because the local pair distribution function g^ (0) approaches one in this limit. 
The gas with infinite repulsion, in turn, has a finite value 

s(oo) = 32 In 2/(3tt 2 ) ~ 0.7491252 (56) 

of the dimensionless contact density, which follows from the asymptotic behavior e(j) = 
e(oo)(l — 41n2/7 + . . .) of the ground state energy (35[. In fact, a similar result also 
applies to bosons in one dimension with an infinite zero range repulsion. This is a Tonks- 
Girardeau gas fl2j , whose equation of state and pair distribution function coincides with 
that of a free Fermi gas with Fermi wave vector kp = ttub- Its momentum distribution, 
however, is very different. At small momenta, it diverges like 1/y/k while for large 
momenta it exhibits a Cs/fc 4 -tail as shown by Olshanii and Dunjko [![. Using equation 
(|25p . which relates the bosonic contact to the pair correlation function at zero distance, 
the asymptotic behavior <7^(0) = 47r 2 /37fj in the limit 7^ > 1 [37| gives rise to a 
universal value 

4 
3"^ 

of the dimensionless contact sb = Cs/kp, which is much smaller than that for infinitely 
repulsive fermions. The finiteness of the dimensionless contact 5(7) in the limit of infinite 
repulsion can be inferred from a simple argument. Indeed, at 7 = 00, the momentum 
distribution n a (k) = f(x) < 1 must be a universal function of the dimensionless variable 
x = k/kp with a power law decay s/x 4 at x 3> 1. The normalization Jdxf(x) = 1 
together with the fact that the power law s/x 4 appears for x = 0(1) due to the absence 
of another momentum scale in the problem beyond kp then implies that s(oo) must be 
of order one. Using equation (1551) . the finiteness of s(oo) immediately implies that the 
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s B (oo) = — ~ 0.1350949 (57) 




Figure 7: Plot of the pair-correlation function (0) times 7, determining the prefactor of the l/g 2 -tail 
of the static structure factor in equation H49B 



energy e(-y) approaches a universal constant e(oo) in the limit of infinite repulsion. As 
will be shown below, a completely analogous line of arguments applies to the 3D repulsive 
Fermi gas with contact interactions. 

Apart from the pair distribution at vanishing distance and the related dimensionless 

contact density s that are shown in figure [__ the homogeneous balanced gas is also 

(2) 

characterized by the combination aiC(ai) ~ 7(7^ (0) that appears both in the pressure 
relation (j33|) and in the asymptotics (l49t of the static structure factor. This combination 
is essentially the interaction energy per particle. For the ground state of the balanced 
homogeneous gas, it can be written in the simple form 

(H') h 2 n 2 

V = ^-7e'(7) (58) 

_____ /"2~) 

where we have used the equations ([_} and (|2"4"|) . A plot of the product 7(7!, (0) is shown 
in figure [__ It exhibits a maximum at 7 = 4.541 . . . with value 1.3266 . . .. Physically, 
this maximum comes about because the interaction energy of the two-component Fermi 
gas becomes small both as 7 <C 1 and for 7 3> 1. In the latter limit, a strong zero 
range repulsion is equivalent to having effectively a Pauli exclusion principle also between 
fermions of opposite spin. The total energy is thus essentially of kinetic origin. The fact 
that e(oo) = 4e(0) shows that the gas with infinite repulsion has the same ground state 
energy as a single-component gas with twice the density, where the associated doubling 
of the Fermi wave vector increases the energy by a factor four. In the limit of infinite 
repulsion, therefore, the balanced gas with total spin zero is energetically degenerate with 
a fully polarized ferromagnetic ground state. Note, however, that despite the fact that 
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both states have the same energy and, in fact, have equations of state p — 2£ that are 
identical at any temperature, their momentum distribution functions arc quite different. 
In the fully polarized state it is the step function of an ideal Fermi gas with an effective 
kp = 2fcp = Tin. The unpolarized state, in turn, is a Luttinger liquid with a tail C/k 4 
that has a universal finite value of the contact density C(oo) ~ 0.749k F . The fact that 
ferromagnetism in one dimension only appears at infinite coupling is consistent with the 
Lieb-Mattis theorem [38] . In the following section, we will show that this result can 
be derived by a simple variational argument, even in the absence of the explicit Bcthc 
Ansatz solution for the ground state energy. 



8. Absence of Ferromagnetism for zero range interactions 

The problem of deriving either necessary or sufficient conditions for the appearance 
of a ferromagnetic ground state of itinerant fermions is a major and - to a large extent - 
unsolved problem in many-body physics. A classic result in this context is the theorem 
of Lieb and Mattis that one never gets ferromagnetism in the ground state of itinerant 



fermions in ID with spin- independent interactions 38, 39]. The Fermi gas with two-body 
contact interactions discussed here is just a special case of this and correspondingly can- 
not have a ferromagnetic ground state. As we have seen above, however, the ground state 
of this system in the infinite repulsion limit becomes degenerate with the fully polarized 
ferromagnet. As was shown by Lieb and Mattis, an infinite coupling constant is indeed 
a necessary condition for the appearance of saturated ferromagnetism in one dimension. 
This result is a simple consequence of a variational upper bound on the ground state 
energy but does not require the explicit solution of the many-body problem by the Bethe 
Ansatz. Remarkably, as will be shown below, a completely analogous argument also ap- 
plies in three dimensions, ruling out saturated ferromagnetism for repulsive fermions at 
finite values of kpa, with a > the scattering length of the pseudopotential interaction. 
This may be viewed as a continuum, low-density analog of a result proven by Kollar, 
Strack and Vollhardt 4(J, 41 1 that there is no saturated ferromagnetism in the half- 



filled fermionic Hubbard model with finite on-site repulsion U unless one adds further 
ferromagnetic exchange interactions, e.g. to nearest neighbors. 

For a balanced, homogeneous ID Fermi gas, the Tan adiabatic theorem in the form 
(|55p together with the fact that the contact density s^j) remains finite at infinite cou- 
pling, quite generally implies that the dimensionless energy e("f) approaches a finite 
value as 7 — > oo. In the absence of the Bcthc Ansatz solution, the Tan relations alone, 
however, give no information about the asymptotic value e(oo). In particular, it is im- 
possible to decide on the basis of the Tan relations alone whether the ground state energy 
e' F ) = 4e(7 = 0) = n 2 /3 of a saturated ferromagnet might be below that of the normal 
Luttinger liquid state with zero polarization beyond a finite critical value of the repulsive 
interaction. A simple variational argument, however, shows that this cannot happen, in 
accordance with the Lieb-Mattis theorem. To see this we note that, for contact interac- 
tions, the fully polarized ferromagnet 

({zz}|* (F) ) =dct<M{ Zl }) |S = N/2) (59) 

is a Slater determinant of the N lowest single-particle states 4n for I = 1, ... TV combined 
with the symmetric state \S = N/2) with maximal spin (note that N = iVf + iVj. is 
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even). Following a well known argument which shows that for just two fermions with 
spin- independent interactions the ground state is always a singlet [38j |. we now choose a 
variational state of the form 



(Ml*™) = Idetft ({zi})| |S = 0) . (60) 

It describes an unpolarized Fermi gas with vanishing magnetization which combines 
the symmetric spatial wave-function formed by taking the absolute value of a Slater 
determinant in ((59)) with the completely antisymmetric spin state \S = 0). Its spatial 
wave-function describes a gas of N bosons with an infinite zero range repulsion. In 
one dimension, this is a Tonks- Girardeau gas, whose energy is known to be equal to 
that of a single-component, non- interacting Fermi gas with wave vector irn (42l [25| . 
The variational state l^var) therefore provides an upper bound e(-f) < tt 2 /3 on the 
dimensionless ground state energy of the unpolarized, repulsive gas. Since l^var) has 
vanishing probability for any of the fermions to be at the same point in space, it can 
be (and actually is) an exact eigenstate only at 7 = 00, where 9^ (0) = 0. The limit 
e(oo) = 7r 2 /3 of degeneracy with a fully polarized ferromagnet can therefore only be 
reached at a point, where the two-component Fermi gas completely fermionizes in the 
sense that there is an effective Pauli principle even for particles with opposite spin. 
Without an explicit solution of the problem, the argument thus rules out saturated 
ferromagnetism in ID at any finite coupling. 

The argument above is instructive but of course it gives no really new information 
because the ground state of ID fermions with contact interactions is known exactly 
by the Bethe Ansatz and ferromagnetism of any kind - not just the simple fully satu- 
rated ferromagnet considered here - is excluded for much more general interactions by 
the Lieb-Mattis theorem. Remarkably, however, the argument can be extended in a 
straightforward manner to the case of fermions in 3D which interact with a repulsive 
pseudopotential. In this case, due to the existence of a two-body bound state for positive 
scattering lengths, the true equilibrium state is an s-wave superfluid that evolves from the 
nontrivial unitary gas at fcpa = 00 toward an ideal Bose-Einstein condensate of dimers in 
the limit kpa — > 0+ [HEEi- The repulsive normal Fermi liquid at positive scattering 



lengths therefore corresponds to a metastable branch of the 3D balanced Fermi gas with 
contact interactions. In the experiment (l9j . the gas is initially prepared in the weakly 
interacting limit kpa <C 1 and then the interactions are quickly ramped to the strongly 
repulsive regime kpa 1 in order to avoid production of bound dimers 43]. An analy- 
sis of the dynamic competition between a putative formation of ferromagnetic domains 
within an RPA-approximation and the formation of pairs indicates that the growth rate 



of the latter instability is always faster (4J). In particular, as shown by these authors, 



anomalies of the kind seen experimentally can be explained in terms of pair formation 
rather than incipient ferromagnetism (see also [IB]]). 

In the following, we will ignore the issue of the dynamics and discuss the question 
whether a Stoner instability might occur in the strongly repulsive Fermi 3D gas with 
contact interactions, assuming a decay to bound pairs happens on a time scale much 
longer than the equilibration time of the repulsive gas (thus, effectively, the state with 
bound pairs is projected out). Since the Tan relations are based on operator identities, 
they also apply to the repulsive branch of the gas. This branch is characterized by a 
dimensionless contact density s r = C r /k F which is determined by the asymptotic decay 
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of the momentum distribution (the subscript in s r refers to repulsive branch). Defining 
a dimcnsionless energy e(k F a) per particle, normalized such that e(0) = 1, via 



AT =5 2^- e{kFa) (61) 

in a manner analogous to equation (|52j) in the ID case, the Tan adiabatic theorem in 
3D @, H reads 

A ' = ^77-^2 -Sr(k F a). (62) 
a{k F a) Z(k F a) z 

From the standard textbook result for the dimensionless energy e(k F a) of the repulsive 



Fermi gas to second order in k F a 17], the contact density in the weak coupling limit 
k F a ^ 1 is given by 

s r (k F a « 1) = ^{k F af + 16(11 1 ~ 2 3 ln(2)) (fcF a ) 3 + ■ • ■ , (63) 

where the leading term just accounts for the mean field contribution. Similar to the 
argument above in ID, the Tan adiabatic theorem (|62p in 3D implies that the energy 
e{k F a) necessarily approaches a finite constant e(oo) in the infinite repulsion limit be- 
cause the dimensionless contact s r of the repulsive gas cannot diverge. In order to prove 
or disprove the appearance of saturated ferromagnetism of the repulsive 3D gas with 
contact interactions, it is necessary to know whether the constant e(oo) that character- 
izes the infinite repulsion limit of the dimensionless energy is above or below the value 
e\ F > = 2 2 / 3 = 1.5874... reached for a fully polarized (and thus ideal!) Fermi gas with 
Fermi wave vector k F — 2 1 ^ 3 k F . The existence of a ferromagnetic instability is suggested 
by a renormalized Hartree-Fock calculation, where a finite magnetization in the ground 
state appears in a discontinuous manner for kpa > 1.054 and full polarization is reached 
already for k F a > 1.112 [l8(. These values are not far from what is obtained by extending 
the perturbative result 

FS(kpa) = ~-k F a - -L (l _ ln(2)) {k F af + ..., (64) 

for the Landau parameter that determines the spin susceptibilit}0 

m */ m (aK\ 

xs - yTf$ (65) 

of a Fermi liquid to k F a values of order one, predicting a divergent Xs with Fq\ c = — 1 
at a finite critical value k F a\ c ~ n/2 + . . .. Moreover, the appearance of saturated 
ferromagnetism in the ground state of the repulsive gas is supported by Monte Carlo 



calculations [21|, |22j which find that the energy of the unpolarized repulsive gas rises 
beyond that of a fully polarized state at a critical coupling k F a\ c of order one. An 
extension of the variational argument used above in the ID case shows, however, that 



1 Recent measurements of the spin susceptibility of the unitary Fermi gas in its normal state have 
been performed by Sommer et al. |4(I . 
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this cannot happen for pure contact interactions. To see this, we choose as a variational 
state for the 3D repulsive Fermi gas with vanishing magnetization a state of the form 
given in . This state has no off-diagonal long range order and is thus orthogonal to the 
superffuid ground state of the 3D pseudopotential at positive scattering length. Its spatial 

1 /3 

wave- function describes a gas of bosons at infinite repulsion n B a = +00. Such a gas is 
known to be effectively fermionized, i.e. its energy per particle is of order H 2 n 2 ^ 3 /m [4^ . 
The tendency to fermionization in 3D is, however, less pronounced in 3D than in ID. 
Specifically, as found by Song and Zhou within a Bogoliubov type Ansatz for the bosonic 
formulation of the problem, the energy of the 3D Bose gas with infinite repulsion stays 
below that of a single-component Fermi gas at the same total density n = ub by about 
20 percent 0|. Thus, there is an upper bound e(oo) < 0.8 on the energy of the 
repulsive, unpolarized Fermi gas which excludes the appearance of a fully saturated 
ferromagnet at finite values of kpa. An important point to note in this context is that the 
argument is restricted to zero range interactions. Moreover, it rules out only saturated 
ferromagnetism for the homogeneous gas not, however, partially polarized states or a 
non-vanishing spin in the ground state of a few-body system. Indeed, the ground state 
of three or five fermions in an isotropic harmonic trap has non-zero angular momentum 
I = 1 [16]. As a result, three fermions with the same spin can be put in the three 
degenerate sublevels of the I — 1 manifold, thus giving rise to a fully polarized ground 
state of three fermions in a trap, as noted by Liu et. al. [49j . The underlying degeneracy 
of the ground state, however, does not extend to larger particle numbers and, indeed, our 
argument above rules out saturated magnetism as the ground state of the homogeneous 
gas. In practice, for ultracold gases, the repulsive branch at positive scattering lengths 
is only metastable. Therefore, a calculation of the energy on this branch alone is not 
sufficient to determine the conditions for ferromagnetism at least as a metastable state 
For example, a variational calculation of the energy of a single flipped spin in the 
background of a fully polarized up-spin Fermi sea with wave vector kp^ indicates that the 
fully polarized state is stable if kp^a > 2.35 [2(|. This calculation, which is analogous to 
that of the spin-flip energy in the Nagaoka state of the lattice Hubbard model discussed 
below, does not account for the finite lifetime of the repulsive polaron. The latter has 
recently been determined by a functional renormalization group approach 50]. It turns 
out that in the 'ferromagnetic' regime kp^a > 1.57, the lifetime h/Y is smaller than about 
bh/ep^. The fully polarized state therefore decays on a microscopic timescale that rules 
out fully saturated ferromagnetism even as a metastable state. 

Within a broader context, the conclusion that a single band model with zero range in- 
teractions shows no saturated ferromagnetism at finite values of the coupling constant is 
in fact not surprising. Indeed, related questions have been studied in considerable detail 
in the context of the 'nearly ferromagnetic' normal phase of Helium 3, where F§ ~ —0.7 
at atmospheric pressure. As discussed by Vollhardt [Hl|, under compression this system 
does not become ferromagnetic with Fq approaching — 1. Instead, it is the effective 
mass m* in equation ([63)1 that diverges, signaling a transition to a solid. Helium 3 is 
therefore an almost localized rather than a nearly ferromagnetic Fermi liquid 51]. Of 
course, the instability to a solid phase requires a hard core part of the interaction, which 
is not present in the case of ultracold atoms. Therefore, it might appear that for pure 
contact interactions a ferromagnetic ground state is not pre-empted by other instabil- 
ities. Approximate calculations of the Fermi liquid parameters for the relevant simple 



22 



Stoner Hamiltonian with an interaction of the form U J n^(x)n^(x) however indicate that 
Fq > —0.63 is bounded from below for arbitrary large values of the repulsion U [H3 |. thus 
excluding a ferromagnetic instability in this simple model. This is consistent with an ear- 
lier argument by Kanamori, which states that the effective repulsion U* that enters the 
Stoner criterion remains finite even in the limit of a diverging bare repulsion U [H3I ] . The 
stability of saturated ferromagnetic ground states has been studied extensively also in 
lattice models of the Hubbard type [54[. In this case, saturated ferromagnetism appears 
for infinite on-site repulsion near half filling, because for a single hole it is energetically 
favorable to move in a background of fully aligned spins 55j . On the basis of variational 
wave-functions it is possible to show, however, that the associated saturated ferromag- 
netism is restricted to the regime near half filling but does not survive at low densities, 
where continuum physics applies. For instance, in a simple cubic lattice the Nagaoka 
state is only stable for hole dopings 8 < 0.24 [56j. A single band Stoner model with 
zero range repulsion is therefore not adequate to account for ferromagnetism, which - in 
practice - requires the presence of flat, degenerate bands and exchange interactions of 
finite range 41 1. 



9. Conclusion 

In this work, we have extended the Tan relations for fermions with contact interactions 
to the situation in one dimension. We have calculated the associated Tan contact in the 
ground state for arbitrary interaction strengths from the exact Bethe-Ansatz results for 
the energy per particle. In particular, we have shown that there is a universal constant in 
the asymptotic decay of the momentum distribution in the limit of very strong repulsion 
which is connected with the saturation of the ground state energy as a function of the 
coupling constant. As in the closely related 3D case, the Tan relations are associated 
with the short range behavior of the one- and two-particle density matrix. As such, 
they may be viewed as complementing the universal features that characterize quantum 
liquids in one dimension in terms of their long-distance, low-energy behavior as Luttinger 
liquids [33|, An important feature of the Tan relations is their applicability to an 
arbitrary state of the system. The assumption of contact interactions, in turn, appears as 



a strong restriction. As emphasized by Braaten [6| and earlier by Zhang and Leggett [57[ , 
however, the Tan relations hold more generally, provided the interaction range rg is much 
smaller than any other characteristic length scale, in particular the average interparticlc 
spacing nT x and the thermal wavelength At- 

A surprising conclusion that appears from our analysis is that a combination of the 
Tan relations with a variational argument rules out saturated ferromagnetism in repulsive 
Fermi gases with contact interactions not only in the - well known - case of one dimension 
but also in three dimensions. In essence, the argument is based on the observation that 
fully saturated ferromagnetism for zero range interactions requires fermionization even 
of particles with opposite spin. This is reached asymptotically for infinite repulsion in 
ID. In the 3D situation, however, the tendency to fermionization is less pronounced. 
Correspondingly, no saturated ferromagnetism appears in this case at finite values of the 
repulsion. Of course, it remains an open problem to study this issue in the presence 
of a finite range interaction, where both the Stoner instability and the transition to a 
solid ground state are possible in 3D [22[ • Moreover, in view of the well known Kohn- 
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Luttinger argument 58j, the ground state of the repulsive gas is likely to be a (p-wave) 
superfluid, whose transition temperature might become of order Tp if kpa ^> 1. 

The possibility to realize strongly interacting Fermi gases with tunable interactions 
in one dimension with ultracold gases in optical lattices 5^, 12 1 opens the chance to 
experimentally test the relations derived here. This is of particular interest since they 
apply even in situations that are still far from the ground state which is hard to reach 
in one dimension due to the typically long equilibration times. 
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Appendix A. Feynman rules 

The matrix elements needed to establish the Operator Product Expansions in sections 
[4] and [6] via matching can be calculated using the Feynman rules for the theory defined 
by the Hamiltonian density in equation ([1]) . They are supplemented by the rules for the 
operator vertices, which can be derived by explicitly Wick-contracting. Alternatively, 
in an approach using standard quantum mechanical scattering theory, one could also 
think of the operator vertices as matrix elements of the operators in the non-interacting 
theory. Note that in this approach, the sums of the diagrams in sections |4] and [6] are the 
analogues to the full matrix element (±p'| £l*_AVt + |±p) of some operator A, where f2± 
are the M0ller operators from scattering theory, such that this a matrix element between 
scattering states that have m-asymptote |±p) and oiti-asymptote |±p'). 

Appendix A.l. Standard Diagrams 

We start with the standard Feynman rules of the theory, needed to evaluate, for 
example, the diagrams leading to the Amplitude A given in (|4"|). 

Momenta. Assign a two- momentum (po,p), where po denotes the energy and p the mo- 
mentum, to all external lines. The two-momenta of internal lines are constrained by 
conservation of energy and momentum at each interaction vertex. 

Propagators. For each internal line of species a with two- momentum (qo,q), assign a 
propagator factor i/(qo — q 2 /2m a + ie). Diagrammatically, the propagation of f-particles 



is symbolized by a solid line (see figure 8(a) ), the propagation of 4,-particles by a dashed 
line (see figure [8(b)] ) . 
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Vertices. For each interaction vertex (see figure 8(c)), assign a factor —ig\. Conserve 
energy and momentum at each such vertex. 



Loop momenta. Integrate over all two- momenta (go, <z), which are not determined by the 
two-momenta of the external lines and two-momentum conservation at the interaction 
vertices, with the measure J dq$dq / (2ir) 2 . 

Appendix A. 2. Operator vertices 

The various operators occurring during the matching contribute their part to the 
matrix elements as well. We list their contributions to the diagrams involving two- 
particle scattering in order of their appearance. In particular, the operator vertices do 
not need to conserve energy or momentum. Let q and — q be the momenta entering, 
where the f-particle carries momentum q and the .J,-particle momentum —q. Let further 
q' and —q' be the momenta leaving, such that the f-particle carries q' and the .(.-particle 
carries —q'. We state explicitly, that the one- particle operators all have in common delta 
functions for the momenta of the particles they don not act on. This delta-function is of 
course no longer there in diagrams only containing one particle, such as those in figure 
[2j Further, for the rules describing one-particle operators of the ^-species, just revert the 
signs of the momenta. 



• ip\{R)^{R + x): 



e i(q-l')R e iqxS(-q + q') 
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Figure A. 8: Basic diagrams for the Feynman rules 
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Appendix A. 3. Example diagram 

Let us evaluate the diagram in figure [3(d)] with incoming momenta ±p and outgoing 
momenta ±p' , which is the part of the one-particle density matrix element that produces 
the non-analyticities, and thus plays a crucial role for the study of the tail of the mo- 
mentum distribution. Using the Feynman rules given above, the diagram is given as the 
integral 



iA(E)iA{E') 



dqpdq 
(2^ 



E-q - 



ieE' - q - 



2m^ 



ie q 



■r 

2m; 



(A.l) 



where we again used the shorthands E = p 2 /2m r , E' = p l2 /2m r . The go-integration 
can be carried out using the Residue theorem. This removes one propagator factor and 
contributes a factor of (— i), leading to 



-iA(E)iA(E') [ ^e l 

J 27T 



= 4m 2 r A(E)A(E') 



1 



1 



E--f- 



2 \p\ 2 ^ \p'\ 



'12 



t|p||x| 

w 



i\p'\\x\ 

w 



(A.2) 



where we have used the Residue theorem again, and carried out the limit e — > 0. For 
x > 0, one has to close the contour through the upper complex plane, for x < one has 
to use the lower one. 
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